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Compact fitting formulas, which contain 
four fitting constants, are presented for 
electron-impact excitation and ionization 
cross sections of atoms and ions. These 
formulas can fit experimental and 
theoretical cross sections remarkably 
well, when resonant structures are 
smoothed out, from threshold to high 
incident electron energies (<10keV), 
beyond which relativistic formulas are 
more appropriate. Examples of fitted 



cross sections for some atoms and ions 
are presented. The basic form of the 
formula is valid for both atoms and 
molecules. 

Key words: electron-impact cross 
sections; excitation; fitting formulas; 
helium; hydrogen; ionization. 

Accepted: November 19, 1992 



1. Introduction 

For applications in plasma modeling and the 
study of energy deposition in matter by charged 
particles, various analytic formulas have been used 
to fit ionization and excitation cross sections, or, of 
neutral atoms, ions and molecules by electron 
impact. The Bethe theory [1,2] provides clear 
guidelines for choosing such analytic formulas at 
high incident energies, E, for electric dipole (El) 
allowed transitions: 

cr{^)==x~\a\nx+b +clx), v/iihx=E/I, (1) 

where / is the ionization potential (or excitation 
energy for discrete excitations), and a , b and c are 
constants characteristic of the target but indepen- 
dent of £, or the scaled incident energy j:. 

The logarithmic term arises from the dipole 
interaction, and hence an El -forbidden transition 
will begin with the b term in Eq. (1). Actually, the 
expansion in Eq. (1) continues after the c term in 



negative powers of x. The Bethe formula is based 
on the first Born approximation. The physical ef- 
fects not represented by the Born approximation — 
such as the distortion of the incident wave, 
polarization of the target charge distribution by the 
incident particle, and the electron exchange effect 
between the incident electron and bound electrons 
in the target— contribute terms between the b and 
c terms in Eq. (1). Hence, expanding the cross sec- 
tion formula beyond those included in Eq. (1) 
becomes a futile exercise unless the additional 
physical effects mentioned above are also included. 
The fitting formulas discussed in this article are 
only valid for El-allowed transitions, since we keep 
the logarithmic term in our formulas. 

According to the Mott scattering formula, which 
includes the electron exchange between the projec- 
tile and target electrons (or an ejected one in our 
case), the interference between the scattered and 
ejected electrons generated by electron-impact 
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ionization leads to a term [2] of the form jc~^ Injc in 
(t{x). On the other hand, the distortion of the 
projectile wave function near the target nucleus 
introduces a term [3] of the form Jc"'xo■(plane- 
wave Born). Moreover, cross sections must vanish 
at the threshold, except for the discrete excitations 
of ions. The terms associated with In x vanish for 
El-forbidden transitions. Various fitting formulas 
[4,5] differ in the choice of functions to represent fl- 
at low X toward the threshold at jc = l, and often 
contain more fitting constants than our formulas. 

To accommodate the leadings: dependence both 
of the electron exchange and the projectile-wave 
distortion mentioned earlier, it is desirable to 
introduce negative powers of x in the fitting for- 
mula. This was achieved by introducing (x+D) in 
the denominator with a fitting constant D, which 
provides additional flexibility. We found that the 
following compact formulas fit known cross sec- 
tions remarkably well throughout the entire range 
of /^£< 10 keV. For a cross section that vanishes 
at the threshold {x = 1), e.g., for the ionization of 
atoms, ions and molecules and for the excitations 
of neutral atoms and molecules, we use 

{r(x)x =S(x) =A \wc + [B \nx + C(^x - l)]/(x +D), 

(2) 

where S(x) is a scaled collision strength and^, B, 
C, and D are fitting constants. Cross sections for 
the El-allowed excitations of ions do not vanish at 
thresholds, i.e., (t(x = 1) > 0. In this case, we use 



ar(x)x=S(x) =A Inx + (B lnx + Cx)/(x +D). 



(3) 



The scaled collision strength, S(x), monotonically 
increases with x in most cases, and the fitting 
constants are better determined using S rather 
than (T itself, particularly when the asymptotic be- 
havior of the cross section at highs; is known either 
through actual measurements or theories such as 
the Born approximation. 

Often, electron-impact cross sections close to 
thresholds are crowded with resonances, mostly 
through the formation of transient multiply-excited 
states or negative ion states. Fitting formulas pre- 
sented in this report are too simple to reproduce 
such resonances. 

In Table 1, we list fitting constants [see Eq. (2)] 
that reproduce electron-impact ionization cross 
sections for some first-row atoms and ions. For the 
hydrogen atom, we fitted to the experimental ion- 
ization cross section measured by Shah, Elliott, and 



Gilbody [6]. For the helium atom, we used the 
experimental ionization cross section by Shah et al. 
[7] and the asymptotic Bethe cross sections by Kim 
and Inokuti [2]. For the ionization of He"^ and Li"^, 
the cross sections measured by Dolder and cowork- 
ers [8,9] and the asymptotic cross sections from 
Ref. [2] were used. Our fitting reproduces the orig- 
inal data within a few percent, except for minor 
local departures mostly near ionization thresholds. 



Table 1. The constanls I,A,B, C, and D [see Eq. (2)] for 
elcclron-impacl ionizalion of atoms and ions. The ionization 
potcnlial / is given in eV, and the fitling constanls /I , B , C, and 
D are in 10''" cm^ 



Atom 



/(eV) 



B 



- C 



D 



H 13.61 0.7576 

He 24.58 0.7326 

He* 54.42 0.06233 

Li+ 75.60 0.08329 



-5.521 5.867 2.948 

-5.117 3.295 2.468 

-0.2982 0.2835 0.327 

-0.4476 0.3225 4.252 



We note that the ionization cross section of the 
hydrogen atom recommended in Ref. [5] is about 
15% higher than the latest experimental value [6] 
at the peak cross section. These two cross sections 
are compared in Fig. 1. 
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Fig. 1. Cross section for the ionization of ihc hydrogen aiom by 
electron impact. Solid circles represent Ihe experiment by Shah, 
Elliou, and Gilbody [6], the dashed curve represents the cross 
section recommended by the Belfast group [5], and the solid 
curve is our fining lo ihe experimental data [6]. 



As an example of El-allowed discrete excita- 
tions, fitting constants [see Eqs. (2) and (3)] for the 
h^^S^lsnp^? («=2-4) excitations of He 
[10,11] and the 2j ^S-*2fi ^P excitation of C^+ by 
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electron impact [12] are listed in Table 2. The He 
excitation cross sections were combined with the 
Bethe cross sections at high incident energies [13]. 
Also, the experimental data on the excitation of 
C^* by Taylor et al. [12] were smoothly joined to 
the plane-wave Bom cross section at high incident 
energies using the shape of the distorted-wave 
Born cross section. The Born cross sections were 
evaluated by the present author using Hartree- 
Fock wave functions for the ground and the 2p ^P 
states as is shown in Fig. 2. Equation (2) should 
also be useful in fitting electron-impact cross 
sections for molecules. 



Table 2. The constants I, A, B, C and D [see Eq. (2)] for the 
Is^ 'S->Unp 'P (n = 1 -4) transition of He and the 2s 'S-»^ ^P excita- 
tion of C^* by electron impact. The excitation energy / is given in eV, 
and the fitting constants^, £, C, and D are in A^ 



Atom Excited state /(eV) A 



B 



D 



He 


ls2p >P 


21.22 


0.3991 


-0.3314 


-0.00325 


-0.012 


He 


U3p ip 


23.09 


0.08978 


-0.0499 


0.01003 


-0.789 


He 


ls4p >P 


23.75 


0.03488 


-0.0286 


0.00571 


-0.343 


C'+ 


2p2p 


8.004 


3.938 


6.373 


0.0061 


-0.999 



2s->2p excitation of C^ 
by electron impact 



o 




f — Fit to exper. 
. Taylor et al. (1977) 
- — Plane-wave Born 
Distorted-wove Born 



where kT is the (plasma) electron temperature in 
eV, and w, and xi are the Laguerre quadrature 
weights and abscissas, respectively. 

Note that the integration formula (8) in Ref. [5] 
has a misprint: ATi there should be replaced by 
(x; +I/kT). Nevertheless, the rates reported in 
Ref. [5] are correct. The fitting constants in Tables 
1 and 2 will lead to cross sections in 10"" cm^ and 
collision rates in cvn?/s. 

For El-forbidden transitions, the logarithmic 
terms in the fitting formulas vanish. The remaining 
two terms with C and D coefficients alone do not 
provide enough flexibility to fit El-forbidden cross 
sections except in unusually simple cases. Instead, 
a power series in (x - 1), 



0fX='ZiAi(x-iy , 



(5) 



2 3 4 5 

ln(E/R) 



where th Ai are fitted coefficients (i^5), should be 
adequate to reproduce most El -forbidden cross 
sections. 

Also, the formulas presented here are not 
suitable for fitting proton-impact ionization cross 
sections at low incident energies, say£~300keV 
or lower, because the analytic form used here 
cannot account for ionization by charge transfer— 
an electron in the target attaching itself to the 
incident proton— which begins to dominate the 
ionization process at E ~300 keV and below. 
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Fig. 2. Collision strength for the 2j ^S->2p ^P excitation of C'+ 
by electron impact. The plane-wave and distorted-wave Born 
cross sections were calculated by (he present author using 
Hartree-Fock wave functions. In the figure, a- is the cross 
section, E is the incident electron energy, ^u is the Bohr radius 
(5.29 nm), and R is the lydberg energy (13.6 eV). 



To calculate the rate coefficients (av) averaged 
over a Maxwellian distribution, we recommend 
that eight-point Laguerre quadrature be used, as 
suggested in Ref. [5]: 

<ov)= 6.692X 10' (kTy'^%(OiS (I +kTx! /I), (4) 
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